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^ : ABSTRACT 

OO ' Context. Past and recent observations have revealed unexpected variations in the far-infrared - millimeter (FIR-mm) dust emissivity 

OA ' in the interstellar medium. In the Herschel spectral range, those are often referred to as a 500 yum emission excess. Several dust 

emission models have been developed to interpret astrophysical data in the FIR-mm domain. However, these are commonly unable to 
' J ' , fully reconcile theoretical predictions with observations. In contrast, the recently revised two level system (TLS) model, based on the 

■ disordered internal structure of amorphous dust grains, seems to provide a promising way of interpreting existing data, 

r ^ ' Aims. The newly available Herschel infrared Galactic (Hi-GAL) data, which covers most of the inner Milky Way, offers a unique 

I opportunity to investigate possible variations in the dust emission properties both with wavelength and environment. The goal of our 

. analysis is to constrain the internal structure of the largest dust grains on Galactic scales, in the framework of the TLS model. 

CLc Methods. By combining the IRIS (Improved Reprocessing of the IRAS Survey) 100 /^m with the Hi-GAL 160, 250, 350, and 500 

fjm data, we model the dust emission spectra in each pixel of the Hi-GAL maps, using both the TLS model and, for comparison, a 
single modified black-body fit. The effect of temperature mixing along the line of sight is investigated to test the robustness of our 
results. 

^ ■ Results. We find a slight decrease in the dust temperature with distance from the Galactic center, confirming previous results. We also 

report the detection of a significant 500 pm emissivity excess in the peripheral regions of the plane (35° < |/| < 70°) of about 13-15% 
of the emissivity, which can reach up to 20% in some HII regions . We present the spatial distributions of the best-fit values for the 
, two main parameters of the TLS model, i.e. the charge correlation length, 1^, used to characterize the disordered charge distribution 

^ ■ (DCD) part of the model, and the amplitude A of the TLS processes with respect to the DCD effect. These distributions illustrate 

' the variations in the dust properties with environment, in particular the plausible existence of an overall gradient with distance to the 

X^"^ I Galactic center. A comparison with previous findings in the solar neighborhood shows that the local value of the excess is less than 

QQ , expected from the Galactic gradient observed here. 

Key words. ISM;dust, extinction - Infrared: ISM - Submillimeter: ISM 

1 . Introduction where Iv(/l) is the brightness, £y is the dust emissivity per hydro- 

gen column. By is the Planck function, Td is the dust tempera- 
The mterstellar medium (ISM) is where matter from diffuse ^^^^^ ^he hydrogen column density. When Equation!!] 

clouds IS brought into stai's. In this context, studying dust evolu- ^^g^ interpret the observed emission, 6, is the mean dust 



O 



- 1—1 

X 



^ tion IS important as variations in the properties of the ISM con- emissivity along the line of sight (LOS). This approximation is 

H , stituents likely affect the star formation process. In addition, un- justified because the medium is optically thin in the far-infrared 

. ^ . derstanding the emission from big grains (BGs) is important in (pjp.) submillimeter (submm). The dust temperature is how- 

interpreting experiments devoted to the observation of the cos- ^^^^ expected to vai-y along the LOS (see Section 4.2.2). The 

mic microwave background (CMB) and its fluctuations. These ^^^j emissivity is usually defined as 
experiments typically operate in the millimeter domain, where 
the cosmological signal is dominated by the Galactic foreground / ^ 

emission, of which dust is one of the components. 

The spectral energy distributions (SEDs) of BGs can be ap- 
proximated by modified black-body emission where ey{Ao) is the emissivity at wavelength Aq, and yS is 

the emissivity spectral index, usually taken to be equal to 

Iy{A) — eyBy{A,Td)NH, (1) 2. In recent years, however, an increasing number of stud- 

ies have found that the emissivity might actually depart from 

* Herschel is an ESA space observatory with science instruments a simple A'^ power-law. Indeed, the emissivity spectral index 

provided by European-led Principal investigator consortia and with im- seems to be temperature-dependent: its va lue has been shown 

portant participation from NASA. to decrease with increasing temperature dPupac et al.L 120031: 
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iDesert et al.L l2008t IVeneziani et aU 1201 Ot iParadis et al.L l201d) . 
This T(j-jS anti-correlation has been extensively debated. Some 
authors claim that this behavior is only the consequence of 
mis-h andling noisy data and/or temperature mixing along the 
LOS dShettv et allilOOOh . However, all the aforementioned stud- 
ies have concluded that the effect was not due to noise, and 
some analyses have included tests that rejec ted temperatu r e mix- 
ing as the cause of the observed trend (iParadis et all l2009t 
iMalinen et al.L 1201 Ih . We note that in some extreme cases, the 
inferred fi values can reach up to 3.5-4 at low apparent temper- 
atures (<10 K), which appears very difficult to reconcile with a 
single value of yS and temperature mixing effects. It appears more 
likely that actual variations in the dust emissivity are at play. 

Emissivity variations also appear to be wavelength- 
dependent with the emission spectrum flattening in the sub- 
millimeter (A > 500 - 600 jum) relative to a single modified 



black-bodv emission (IWrieht et al.L 199 


I Reach et all 1995 


Finkbeiner et al.'. 


1999^, 'Gallianoetali 


2005; iPaladini et al. 



2007; Paradis et al., 2009, 201 Ij). They als o app ear to be 
environment-dependent (see, for instance, 0> l2005 u for a re- 
view), with an absolute emissivity value in the FIR {A < 
500 - 600 nm) that is higher in colder and denser environments 
('Stepnik et al.L l2003t IParadis et al.L l2009t iPlanck collaborationi 
l2011u). Moreover, the analysis of the Planck data in the solar 
neighborhood and both the Larg e and Small Magellanic Clouds 
(iPlanck collaborationi |201 lellm|) indicates that the flattening of 
the emissivity spectrum could be related to the metallicity, be- 
cause it is more significant in low metallicity environments. 
From the astrophysical point of view, grain aggregation in the 
core of dense clouds could, in principle, explain the observed 
absolute FIR emissivity excess, if one assumes that aggregates 



emit more FIR radiation than individual grains. Nevertheless, 
this interpretation is only partially satisfactory, since grain ag- 
gregation may affect only the absolute values and not necessarily 
the shape of the emissivity, and the aforementioned excess has 
only been reported for some Galactic molecular clouds, while in 
the extragalactic context, such as in the Large Magellanic Cloud 
(Paradis et al.t,,201 lb) , it has not been found at all. 



In this work, we study the FlR/submm data obtained for 
the first and fourth Galactic quadra nts through the Hersc hel 
infrared GALactic (Hi-GAL) survey (Molinari etall l2010allbl) . 
and demonstrate the existence of a submm emissivity excess. 
Thanks to the sub-arcmin resolution, comprehensive sky cover- 
age, and broad wavelength range spanned by the Hi-GAL data, 
we can for the first time perform a systematic investigation of 
such an excess in our Galaxy. In addition, the Galactic plane is 
ideal for conducting such an analysis because it contains a va- 
riety of different environments, from diffuse clouds to compact 
sources, and from cold to warm structures, which can potentially 
harbor different dust properties. We also show that we are able 
to reconcile the observed excess with theory predictions using 
the two level system (TLS) model, which provides a sophisti- 
cated analytical description of the disordered internal structure 
of amorphous dust grains. 



The paper is organized as follows. In Section|2] we introduce 
the observational data that we use in this work. Sections[3]and|4] 
focus on the variations in the emissivity spectral shape and the 
origin of these variations. Section|5]provides our conclusions. 
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Fig. 2. From left to right: 500 fim emissivity excess (in percent) and dust temperature derived from a modified black-body emission, 
and the results of the TLS modeling for dust temperature (Ttls) in intensity of the TLS processes (A), correlation length (Ic) in 
nm, and reduced 
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Fig. 3. 500 jum emissivity excess (panel A), dust temperature (panel B), A parameter (panel C), and Ic parameter (panel D), as a 
function of Galactic longitudes. The dashed line shows the solar neighborhood value when normalized to ca j=475. The vertical 
error bars correspond to the 1-cr standard deviation. The blu e arrows show the tangent directions of the spiral arms at l=-50°, -33°, 
-21°, and 31° (,Vallee>,2008,) and of the 3 kpc ring at |/|=23° (iDame & Thaddeusi l2008h . 



2. Data 

The reference database for our investigation consists of the Hi- 
GAL dataset, a collection of 56 tiles of 2°x 2°, with the longi- 
tude range -70° < 1 < 70°, and an extension of 2° in latitude, 
following the Galactic warp. The Hi-GAL PAC S (70 and 160 
;<m) and SPIRE (250 350 , and 500 fim) dat a dPilbratt et al.L 
120 1 Ot iPoglitsch et all I20T0I: iGriffin et al.Ll2010h were processed 
using the ROMAGAL software, described in iTraficante et al.l 
(l2011h . We also used ancillary IRIS (Improved Reprocessi ng of 
the IRAS Survey, see lMiville-Deschenes & Lagacheil20()5h data 
to constrain the peak of the BG emission. We corrected the gains 
and zero-level in the SPIRE and PACS images by applying gains 
and offsets derived from the comparison to a combination of the 
Planck-HFI (High Frequency Instrume nt) and the IRIS obse rva- 
tions, using the same procedure as in [Bernard et al.l (l2010h . A 
forthcoming paper will be dedicated to the inte rcalibration of 
the H i-GAL data with existing photometric data ([Bernard et al.L 
201 lb. A ful l descr iption of the Planck mission i s provid ed in 



Tauber et al.1 ( 1201 Ol) and the iPlanck collaboration! (1201 lab . The 
Hi-GAL data (0 = 6" to 37" from 160 to 500 yum) were con- 
volved to a 4' angular resolution to match the resolution of 
the IRIS data. In addition, each tile was re-projected using the 
HEALPIx pixelization scheme (Hierarchical Equal Area isoLat- 
itude Pixelization)' with nside=2048, which corresponds to a 
pixel size of 1.7', i.e. adequate for Shannon sampling of the 
adopted resolution. We used the surface area of the intersection 
between the HEALPix and the native pixel as a weight to re- 
project the data. After reprojecting, all the fields were combined 
into a single HEALPix file covering the entire inner Galactic 



plane. We adopted the same absolute calibration accuracy as for 
the SDP fields in the PACS data, i. e. 20%, and 13.5% for IRIS 
100 fim. The calibration of the SPIRE data was recently revised 
leading to an accuracy of 7%. As a conservative approach, we 
did not include the PACS 70 i^m data in the analysis because 
at this wavelength we are probably most sensitive to emission 
from very small grains ( VSGs), which are stoc hastically heated 
by photon absorption ( Compiegne et al.L I2OIOI) . This contribu- 
tion can only be roughly evaluated and can severely bias the 
derivation of the BG temperature (especially for Td < 25 K). 



3. Emissivity excess at 500 //m 

A break in the observed emissivity law around 500 yum has 



been identified in the outer Galaxy (iParadis et al 



the an alysis of DIRBE (Hauser, 1993), Archeops ( Benoit et al.L 



2009) from 



http://lambda.gsfc.nasa.gov/ 



l2002h . and WMAP (.Bennett et al.i, aOOl) data, as well as in the 
Large Magellanic Cloud bas ed on Herschel data dGordon et al.L 
1201 Ot iGaHiano et all |201 ih . However, in the two Hi-GAL 
Science De monstration Phase (S DP) fields (centered at 1=30° 
and 1=59°), iParadis et alJ (l2010h did not find any significant 
emissivity excess at 500 fim, even when taking calibration un- 
certainties into account. The goal of the present work is to in- 
vestigate the existence of this break. For this purpose, we fit the 
Hi-GAL/IRIS emission in the range 100 fim < A < 350 yum in 
each pixel with a single modified black body, adopting a least 
squares fit method. We assume that calibration uncertainties are 
distributed following a Gaussian function. The dust temperature 
is set as a free parameter We then compare the observed emis- 
sivity at 500 yum (e™) with the value predicted by extrapolating 



Paradis et al.: Detection and characterization of a 500 yum excess in the Galactic plane 



5 




10 20 30 40 



Ta (K) 

Fig. 4. Temperature derived using the TLS model (Jtls), as a 
function of temperature derived from a modified black-body fit 
with;6=2(r/;,). 



the best-fit emissivity power-law to 500 jjm. (e^ol))- We define the 
500 ixm excess as the quantity 



"500 

To perform the fit described above, we need to adopt a value for 
the emi ssivity spectral index fi. A recent analysis of the Planck 
data tPlanck collaboration! l201 1^ has found, for 100 fim < A < 
850 um, a m e dian fi value of 1 .8 over the entire sky. In contrast, 
iParadis et ai] (1201 Cft obtained, for 100 /zm < i < 500 yum, a 
median yS value of 2.3 in both of the Hi-GAL SDP fields. In the 
light of these results, we decide to keep /? fixed and equal to 2. 
We note that this assumption is a matter of how the excess is 
defined and only affects the amplitude of the excess and not its 
spatial distribution. 

In Figure [1] (panel A), we show the pixel distribution of the 
emissivity excess at 500 ;um (esoo) derived from the fitting pro- 
cedure. The distribution is clearly not centered at zero (as would 
be the case for pure noise), but instead visibly skewed towards 
positive values. The median of the distribution is 0.8%, with a cr 
of 4.2%. Excess values larger than 13.4% exceed the 3-cr stan- 
dard deviation of the distribution and are significantly larger than 
the calibration uncertainty of the SPIRE data (cf. Section|2|i. 

We emphasize that this result, i.e. the existence of a signifi- 
cant excess of emission at 500 ^m, does not depend on the par- 
ticular value adopted for fi in the modified black-body fit. If we 
repeated the analysis with a different spectral emissivity value, 
for example (3 equal to 1.8 instead of 2, the shape of the excess 
distribution would remain unchanged, although it would shift 
towards more negative values, suggesting that, for the majority 
of the pixels, the model overestimates the data. Incidentally, we 
note that this also corroborates the hypothesis that, at least along 
the inner Galactic plane, /? is indeed larger than 1 .8. 

Figure |2]illustrates the spatial distribution of the 500 fim ex- 
cess across the Galactic plane. Inspection of the map reveals a 
more pronounced excess in the peripheral regions ((35° < |Z| < 
70°)) than in the central ones (|/| < 35°). This effect also appears 
in Figure [3] panel A, where we have plotted the median value 
of the excess at each Galactic longitude, and its associated 1-cr 
standard deviation. In particular, for -70° < I < -60° the excess 



becomes significant, contributing to as much as 16 - 20% of the 
total emissivity, and even reaching values larger than 23% to- 
ward some HII regions. Regions in the proximity of the Galactic 
center, where a ring of dense a nd cold clouds has b een identi- 
fied at the Herschel resolution jMolinari et al.L 1201 lb . however, 
do not exhibit any excess. The spatial distribution of the excess 
does not display any particular symmetry either: for instance, 
the significant increase in the fourth Galactic quadrant around 
l<-50° is not mirrored by a corresponding increase within the 
same longitude range in the first quadrant. Remarkably, the mor- 
phological behavior that characterizes the 500 yum excess is - 
per se - an indication that this effect is not related to calibration 
uncertanties, as these would have an equal impact on different 
regions of the plane. In the specific fields of 1=30° and 1=59°, 
the median value of the excess is 1% and 3%, which does not 
correspon d to any signi f icant detection, in agreement with the 
findings in Paradis et a 

The submm excess could, in principle, be induced by the fit- 
ting procedure. To test this hypothesis, we carried out two tests: 
first, we included in the fit all data points for 100 < A < 500 
jum; second, we performed the fit by replacing the 350 //m data 
point with the 500 //m data. Both tests provided larger reduced 
than the original fit (on average, by a factor of 2.3 and 2.6), 
thus supporting the idea that this excess has an astrophysical ori- 
gin. 

In the following, we investigate in detail the 500 excess 
in the framework of the TLS model, by analyzing variations in 
the model parameters. 

4. Towards an understanding of the nature of the 
excess 

4.1. The TLS model 

We now prov ide a brief description of the TLS model. We refer 
the reader to Menv et al. (2007) for a fuller description. 

Our revised TLS model is based on both solid state physics 
and laboratory measurements, and provides an accurate descrip- 
tion of the physical properties of amorphous grains. It consists 
of two parts: (1) the disordered charge distribution (DCD) one, 
which describes the interaction between the electromagnetic 
wave and acoustic o scillations in the disordered charge of the 
amorphous material (IVinogradovi 1 1 96^ ISchlomannl 1 1 964 ): and 
(2) the actual TLS part, which takes the interaction of the elec- 
tromagnetic wave with the simple distribution of an asymmetric 
double-well potential into account. This second part corresp onds 
to a theory orig i nally developped by Phillies ( 19721 Il987h and 
lAnderson et alj (|1972|) . The first effect (DCD) is temperature- 
independent, occurs on the grain scale and is the most dominant 
effect in the FIR. Moreover, it has two types of asymptotic be- 
haviors at both lower and higher frequencies than wq: fy A^^ 
in the short wavelength range, i. e. oj > wq, and ey oc A^'^ in the 
long wavelength range, i. e co < wq. This DCD effect is charac- 
terized by the so-called correlation length Ic, which determines 
the wavelength at which the inflection point between the two 
asymptotic behaviors occurs. This parameter is defined in the loq 
term by 

1^0 = vt/lc, (4) 

where Vt is the transverse sound velocity in the material. The 
TLS effects, which consist of three mechanisms (a resonant ab- 
sorption and two relaxation processes, thermally activated), are 
instead temperature-dependent, take place on the atomic scale, 
and start to be important in the submm domain, becoming the 
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Fig. 5. Average emissivity spectra in five bins of temperature and 
their l-cr uncertainty: 16 K in red (solid), 20 K in dark blue 
(dash), 24 K in light blue (dot), 28 K in green (dash-dot), and 
32 K in orange (dash-dot-dot). The emissivities have been com- 
puted using the best-fit with the TLS model (top panel) and with 
a modified black body with j8=2 (bottom panel). An emissivity 
power-law in A^^ is overplotted as a continuous gray line, for 
comparison. 



predominant effects in the mm wavelength range. The amplitude 
of the TLS effects with respect to the DCD part is determined 
by the intensity parameter denoted with A. This means that all 
TLS processes are multiplied by the same intensity A. The total 
emission deduced from the TLS model (/,„,) is then the sum of 
the emissions coming from the DCD (Idcd) and TLS processes 
(Itls ), such as 



Itot - Idcd + ^ Itls ■ 



(5) 



iParadis et alJ (l201 lah performed a comparison between the 
model predictions and astrophysical data to determine the model 
parameters. They analyzed both the FIR/mm SED of the solar 
neighborhood using FIRAS and WMAP data, and of Galactic 
compact sources observed with the Archeops experiment. The 
purpose of their analysis was to constrain the four key parame- 
ters of the TLS model, namely: (1) the charge correlation length 
Ic for the DCD absorption; (2) the amplitude A of the TLS ef- 
fects, with respect to the DCD process; (3) a parameter ca de- 
scribing the double-well potential for one of the TLS processes; 
and (4) the dust temperature Ttls ■ The combined fit of the SEDs 
for the solar neighborhood and the compact sources provided a 
general description of Galactic dust in terms of these parameters. 



whose best-fit values are /c_j=13.4 nm, A s=5.8, and ca,.s =475 (see 
also Table [U, which were referred to as the standard values. 

In the following, we repeat the fitting procedure described in 
Section [3] using this time the TLS model rather than a modified 
black body. The best-fit values for the TLS parameters derived 
from the fit are then compared to the standard ones, and those ob- 
tained for the solar neighborhood and compact sources. Before 
describing the results we obtained, a few considerations need to 
be made. First, we note that the ca parameter is, to first order, de- 
generate with the A parameter In addition, ca cannot be tightly 
constrained in our study because we are limited to wavelengths 
shorter than 500 jjm. For this reason, following IParadis et al.i 
(1201 lah . we set ca to be equal to the standard value, (ca,s = 475). 

As for Ic, we adopt 36 nm as an upper limit. For 1^ > 36 nm, 
the DCD process indeed reaches an asymptotic behavior in A^^ 
from FIR to mm wavelengths. Values of Ic of the order of a few 
nanometers, combined with either small or null values of A, al- 
low p to be instead close to 4. The slope of the spectra increases 
with decreasin g l^, in the absence o f significant TLS processes. 
As analyzed in lParadis et al.l (1201 lal) . for the range 100-550 fim, 
A=0.1 and lc-1 nm (30 nm), the parameter /3 is expected to be 
constant with a value of 3.4 (2.0), regardless of the dust temper- 
ature (see their Figure 7, middle panel). However, for A=10 and 
nm (30 nm), /3 varies from 2.45 (^1.85) to 1.4 (1.35) for 
temperatures between 10 K and 55 K. 

4.2. Dust properties along the Galactic plane 
4.2.1. Method 

To minimize the computing time, we pre-calculated the bright- 
ness in the IRIS 100 jum and Herschel bands using the TLS 
model, taking into account the color-correction to be applied to 
each instrument. In this way, we generated a multi-dimensional 
grid, for different Ttls between 10 K and 40 K (with a 0.5 K 
step), le in the range 1-36 nm (with a 1 nm step), and 50 values 
of A between and 50. We then compared the brightness ex- 
pected from the TLS model with the observed brightness in each 
pixel of the map. The;^'^ value was computed for each value of 
the grid, and we chose the value of the parameters (p) that min- 
imizes the;^'^. Calibration uncertainties were included in the fit 
as described in Section [3] To allow interpolation between indi- 
vidual entries of the table, the best-fit parameter values (/?*) are 
computed for the ten smallest values ofx^ using 



(6) 



Since Planck data are not yet available, we are unable to place 
tight constraints on the A parameter using only Herschel data. 

A normalization of the dust emission predicted by the TLS 
model is required. We adopt the normalization 



. _ Itls x < hbs > 

^TLS ,norm — , 5 

< iTLS > 



(7) 



where Itls and Itls, norm are the emission spectrum from the TLS 
model before and after the normalization, while < I„hs > and 
< Itls > are the emission averaged in wavelength between 100 
jum and 500 yum, corresponding, respectively, to the observations 
and the TLS model. 
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Fig. 6. Left panel: a paramete r from the iDale et alJ (l200ll) model. Right panel: reduced deduced from the comparison between 
data and the lDale et all (1200 Ih model. 



8 



Paradis et al.: Detection and characterization of a 500 yum excess in the Galactic plane 



4.2.2. Dust temperature and emissivity 

The spatial distribution of the dust temperature, as well as the 
reduced obtained through the TLS modeling are presented 
in Figure |2l The temperature histogram is shown in Figure [1] 
(panel B). We note that the reduced values are systematically 
smaller than 1, suggesting that the uncertainties in the observa- 
tional data are probably overestimated, particularly for the PACS 
data. The median temperature along the Galactic plane is 18.6 K, 
while the Galactic ce nter exhibits temperatu res in the range 18- 
25 K. As reported in iMolinari et al.l (1201 Ih . some cold clumps 
can be found in this region, with temperatu r es of about 16-18 
K. However, with respect to iMolinari et al.l (1201 ih . we do not 
find any evidence of temperatures above 25 K. This may be for 
several reasons. First, our temperature map was generated at 4 ' 
angular resolution, instead of 25" as in iMoUnariet all (1201 Ih . 
Secondly, the 70 fim was not included in our analysis, which 
might have caused the underestimation of the warmest dust tem- 
peratures. 

In general, temperatures appear to be higher in the cen- 
tral regions of the Galaxy (Ttls^19A K for |/| < 35°), 
and to decrease with increasing longitude (r7-£s=18.0 K for 
35° < |/| < 70°), as one can clearly see in Figure [3] (panel 
B). This general beha vior has been seen before u sing DIRBE 
(|Sodroski etar.T99?), FIRAS ( Reach et all [19951) . and Planck 
data (Planck collaboration, 201 1^. Temperatures span the range 
16 K - 26 K (Figure [J] panel B). A the IRIS resolution, the 
Galactic plane remarkably does not show dust temperatures 
below 15 K. Some warm spots with Ttls > 30 K are, however, 
visible in the peripheral regions, probably owing to HII regions. 
These warm spots also show the largest 500 /im emissivity 
excess. The comparison between temperatures estimated with 
the TLS model and those derived from the modified black-body 
fit in the range 100-350 ^um with P=2 (Tfn) is provided in 
Figure |4] As described in ' Paradis et al.l ( 1201 la l). there is good 
agreement between the two methods up to ^25 K, where a 
departure from a one-to-one correlation becomes noticeable. 
The dust temperature is only taken into account in the Planck 
function in the case of a modified black-body fit, whereas the 
TLS processes are temperature-dependent and increase with 
temperature. As a consequence, the difference in the temperature 
contribution in the two methods could explain the departure in 
the TxLS-Tfit correlation at high temperatures, and especially the 
faster increase in Tjls compared to Tfn for temperatures higher 
than 25 K. The outlier points (T/,-, > 20 K) correspond either to 
pixels with a low quality for the modified black-body fit, or 
to edges of the map. 

In Figure |5] we have plotted the average dust emissivities, 
for 100 fim < A < 500 yum, in five dust temperature bins cen- 
tered at 16, 20, 24, 28, and 32 K. The plots have been created 
for both the TLS model and a modified black-body fit with j6=2. 
In the first case (TLS model), we computed the dust emissivity 
spectra using the best-fit model parameters (Ttls, A, and 1^) for 
each pixel, that we averaged per temperature bin. We derived 
the averaged emissivity in the photometric bands by applying 
the adequate color-correction for each instrument. In the second 
case (modified black-body fit), we divided the brightness com- 
ing from the data by the Planck function at the temperature Ty,,. 
The emissivity spectra were normalized to unity at 100 jum. The 
uncertainties in the emissivity were computed using 

Ae„ _ ABATd) My 



The relative error ABy(Td)/By(Td) is related to the error on 
(ATd) as 

Afiv(rrf) ^ w-/^ ATd 
Bv(Td) - I Td 

The figure highlights the existence, for both models, of a 
clear trend in the emissivity spectra with the dust temperature, 
with the spectra being flatter with increasing dust temperature. 
This observed flattening behavior with Td cannot be accounted 
for by the calibration error at 500 fim. If the brightness at 500 
fim were systematically off by 7% owing to calibration uncer- 
tainties, all emissivity spectra would be affected but this should 
not induce a change in the dust temperature. Another effect has 
to be considered. Given that we conducted our analysis along the 
Galactic plane, where temperature mixing is expected along the 
LOS, one may argue that we can introduce an artificial bias by 
fitting the data points with a single dust temperature. To address 
this question, following Dale et al. (2001), we computed the pre- 
dicted brightness in each IRIS and Herschel band for different 
values of the interstellar radiation field intensity, Xjsrf, relative 
to the intensity of th e ultraviolet field in the solar neighborhood 
(Mathis etal.V1983). 

The Dale et al. (200T|) model describes a temperature distri- 
bution using the concept of local SED combination, and assumes 
a power-law distribution of dust mass subjected to a given heat- 
ing intensity dMdiXjsRp) 

dMdiXjsRF) XJsupdXjsRF, (10) 

where a and X/srf are in the range l< a < 2.5 and 0.3< Xisrf < 
10^. With this representation, the diffuse medium has a close to 
2.5, whereas active star-forming regions have a of 1. We esti- 
mated the emission spectra {I"!°'^{Xisr f)), for different values of 
a, using the DustEM package (Compieg ne et al.L l2008l 1201 Ih 
and assuming [3-2. We then sum med these contrib utions over 
the same X]srf range proposed bv lDale et al.l (1200 ih 

jtot _ ^ j=\ ^ y^ISRFj) X A^^^p . 

" ~ y88 y50 y-aj " ^^^^ 

Zj/=1 Zjj=\ -^ISRFJ 

We then compared the observational data with the predicted 
SEDs, and searched for the optimal a value that minimizes the 
difference between the two. As a result of this process, each pixel 
of the map has an emission that incorporates temperature mix- 
ing along the LOS. Figure |6] is an illustration of the pixel-to- 
pixel map of the recovered a values. The central regions of the 
Galactic plane are characterized by values of a in the range [1.9, 
2.5], while for \l\ > 30°, a appears to be systematically equal to 
the asymptotic value of 2.5. At the same type, inspection of the 
reduced spatial distribution (also shown in Figure |6]l reveals 
that the SEDs for most of the peripheral regions are not well re- 
produced by the lDale et al.l (2001) model. In these regions, tem- 
perature mixing is probably less important than in the central 
regions of the Galaxy. 

We considered only pixels with reduced < 1.5. Modeling 
with a reduced > 1.5 does not give a satisfactory fit. The spa- 
tial distribution of the reduced is provided in Figure |6] This 
selection removes 39% of the pixels, and in particular the pe- 
ripheral parts or the map, corresponding essentially to the cold- 
est parts of our map. We then proceeded as before, with the 
di fference tha t we tr eated the brightness derived from the best- 
fit iDafeetin ( 1200 Ih model for the selected pixel (61% of the 
Galactic plane) as "real" data: we performed ax^ minimization. 
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100 1000 

X (|im) 

Fig. 7. Emissivity spectra in five bins of temperature and their 
1-cr uncertainty, assuming temperature mixing along the LOS 
(with the use of the Dale et al. (2001) model): 19.1 K in red 
(solid), 22.3 K in dark blue (dash), 25.5 K in light blue (dot), 
28.7 K in green (dash-dot), and 31.9 K in orange (dash-dot-dot). 
The emissivities have been computed for the best-fit TLS model 
(top panel) and a modified black body with y6=2 (bottom panel). 
An emissivity power-law in A'^ is overplotted for comparison in 
continuous gray line. 



using the TLS model and a modified black-body emission model 
with y6=2, to derive the best-fit parameters and the emissivity 
spectra for different temperature bins. Since the coldest regions 
had been rejected from the analysis, the temperature bins ranged 
from 19.1 K to 31.9 K. 

If temperature mixing did indeed account for the flattening 
of the spectra with increasing dust temperature, we would expect 
to observe the same type of behavior with this test. However, 
as shown in Figure |7] the spectra do not show any sign of a 
dependence on dust temperature, regardless of the model used 
to derive the dust emissivities. Thus, temperature mixing is not 
likely to be responsible, and we are confident that the emissiv- 
ity changes are instead related to the intrinsic properties of dust 
grains. 

4.2.3. Model parameters: A and Ic 

Unlike in iParadis et al.l (l201 lah . we explored both cold and 
warm/hot regions of the Galaxy, to systematically investigate the 
variation in the TLS model parameters with environment. 

Figure[8]shows the median values derived from this study for 
the Ic and A parameters in the case of, respectively, the whole 



Galactic plane region covered by Hi-GAL and, the central and 
peripheral parts of the plane. We have also displayed, for com- 
parison, the best-fit values deduced by Paradis et al. (201 la) for 
the solar neighborhood, for Archeops compact sources, and for 
the combination of both. In this analysis, the uncertainties in the 
parameters are assumed to be equal to their standard deviation. 
Since the error bars in the data are likely to be overestimated, 
as already discussed in Section 14.2.21 the parameters cannot be 
strongly constrained. 

The distribution of points in the figure reveals that differ- 
ent environments are likely associated with different dust grain 
properties. From the fitting procedure discussed in Section 4.2.1, 
we estimated that the most representative values of Ic and A for 
the Galactic plane are 7.5 nm and II.O. For Ic, this is close to 
the value derived for Archeops sources, but smaller than for the 
solar neighborhood. We also note (Figure [1] panel D) that Ic is 
globally characterized by values much lower than the standard 
one {Ic^s - 13.4 nm). Both of these results are indicative of a less 
organized charge distribution for dust grains in the plane than 
for dust particles in the solar neighborhood. As for the ampli- 
tude of the TLS processes. A, this is larger by a factor of 1 .9 in 
the Galactic plane than the standard value (A, = 5.8). We note, 
however, that in the solar neighborhood, probed by FIRAS ob- 
servations, and in Archeops compact sources, ca differs signifi- 
cantly from the standard value. Therefore, to be able to perform 
a proper comparison between these results, we needed to nor- 
malize A to ca.j=475. Accordingly, we obtained A equal to 4.8 
and 10.9 for, respectively, the solar neighborhood and the com- 
pact sources. Again, A in the Galactic plane is different from the 
value derived in the solar neighborhood. This result indicates 
that grains within the plane might have a larger degree of amor- 
phization than in the solar neighborhood on the atomic scale, 
where the TLS processes occur 

The spatial variation across the Galactic plane of Ic and A is 
provided in Figure|2] Important variations in the two parameters, 
on both small scales (pixel to pixel) and large scales (from the 
central to the peripheral parts of the Galaxy), are visible. The 
A parameter exhibits high values in regions with significantly 
high 6500 (especially in peripheral regions, see Figure [3j panel 
C), which suggests that the excess does not depend on Ic, but 
rather on A. This explains the large fluctuations in the Ic pa- 
rameter in the peripheral regions (see Figure |2]l, which induce 
large standard deviations in the longitude profile of this parame- 
ter (see Figure[3]). In contrast, Ic is more accurately characterized 
in regions of negative excess. In that latter case, variations in the 
spectral shape of the SEDs can only be accounted for by vari- 
ations in the correlation length, since this parameter allows us 
to produce emission spectra with a spectral index larger than 2. 
This explains why the Ic standard deviation is smaller closer to 
the Galactic center The correlation between the 500 fj.m excess 
and the A parameter is shown in Figure |9] Disorder in the grain 
structure on the atomic scale, i. e. the TLS part of the model, thus 
appears to be the main contributor to the emissivity variations at 
500 jum. The Ic parameter is indeed almost constant along the 
Galactic plane, as one can discern in Figur^D (panel D). 

The A values (or excess values) are higher over most of the 
inner Galactic plane than in the solar neighborhood. If the origin 
of the variations is a general Galactic gradient, this would in- 
dicate that the solar neighborhood is an atypical place. It might 
also suggest that the spatial distribution of the excess is not de- 
scribed by simply an overall Galactic gradient but a more com- 
plicated structure. 

The relation between A and Ic is illustrated in Figure[TOl The 
plots are very noisy, especially for the whole and the peripheral 
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Fig. 8. Best-fit model parameters Ic and A and their 0.5-cr 
standard deviations, derived from the TLS model for different 
Galactic environments. The arrows indicate values of the param- 
eter to first order, when fixing ca to the standard value equal to 
475. 




Fig. 9. A parameter derived from the TLS model as a function 
of the 500 i-im emissivity excess. The overplotted contours, with 
levels at 50, 200, 350, 600, and 1000, represent the density of 
pixels. 



regions of the Galactic plane. Nevertheless, they reveal a weak 
anti-correlation between the two parameters, which disappears 
at large values of A. We also see that, in the central regions of 
the Galactic plane, most of the points cluster in the region de- 
limited by A ^ 6 - 10, whereas the peripheral regions is rather 
accurately characterized by A ^ 10 - 22. This shift in A is likely 
associated with variations in the dust properties, while the disap- 
pearance of the anti-correlation between A and Ic confirms that 
spectral variations at wavelengths equal to or larger than 500 yum 
could be entirely due to variations in the intensity of the TLS 
process. A significant flattening of the spectra in the submm/mm 
domain is indeed expected if the TLS effects are responsible for 
the observed excess. We remind the reader, however, that longer 
wavelength observations are required to confirm the results ob- 
tained for this parameter We note that the calibration uncertainty 
in the 500 /urn measurement does not affect the spatial distribu- 
tion of the TLS parameters. The calibration error would indeed 
only induce a variation in the absolute value of A. For instance, 
a decrease in the SPIRE 500 jim data by 7% would cause a cor- 
responding decrease in A by a factor of 1.6, but the peripheral 
regions of the plane would still have the largest values of A. As 
for Ic, its median value - and standard deviation - would be lower 
(Ic -5.8 nm) than the regular median value (Ic -7.5 nm), even 
if the change in the dust temperature is as small as 0. 1 K. The 
difference of 1 .4 K between the average dust temperature in the 
central and peripheral regions of the Galaxy would be preserved 
(see Section l4.2.2b . 

4.2.4. Relationship witli tine Galactic spiral arms 

According to lGreenberg & Lil (Il995h . the spiral arms regions and 
the interarm regions might have different dust properties. This 
could explain the variat ions with Galactic longitude observed in 
the extinction data of iGao et alJ (l2009l) . The latter authors ar- 
gued that these variations could be caused by the larger grains 
produced by coagulation, a higher dust d ensity and a s t ronger 
radiation field in the spiral arms regions. Bernard et al.l (1201 Ol) 
also showed the existence of warmer dust in regions of the LOS 
corresponding to the intersection with spiral arms. In Figure [3] 
we have marked the tangent positions of the spiral arms taken 



from IValled ( l2008l) (l=-50°, -33°, -21°, and 31° for the Crux- 
Scutum, Norma, Perseus, and Scutum-Crux arms, respectively), 
as well as the position o f the 3-kpc molecular ring defined in 
iDame & ThaddeusI (l2008h (|;| = 23°). 

From Figure |3] one can observe a direct link between the 
tangent positions of the spiral arms with some peaks in the dust 
temperature profile. Variations in the 500 /im emissivity excess 
linked to the arms are also visible but are less significant than 
the general trend with distance from the Galactic center For the 
A parameter, the relationship with spiral arms is less obvious. A 
hint of a relationship can be discerned between the spiral arms 
and the correlation length, which is enhanced in these regions, 
although by an amount that is smaller than the standard devia- 
tion. 

Variations in the optical properties predicted by the TLS 
model replace or are complementary to the predicted variations 
in emissivity caused by either aggregation or fluffiness. Dust 
growth is expected in dense molecular clouds, and aggregation 
likely produces both an increase in emissivity with respect to the 
gas column (for a given dust to gas ratio) and a decrease in dust 
temperature. The increase in the absolute emissivity with respect 
to the gas cannot be measured here because we do not directly 
compare the emissivity with the gas column density. From dipole 
discrete approximations calculations, the emissivity increase for 
aggregates is expected to be roughly constant in the range 100- 
500 yum. This does not change the spectral shape of the emis- 
sivity and cannot reproduce any significant 500 fim emissivity 
excess. Similarly, the changes in apparent dust temperatures are 
difficult to interpret in terms of aggregation at the scales ana- 
lyzed here, and are probably determined mainly by variations in 
the local radiation field intensity, which may obscure any tem- 
perature decrease caused by dust aggregates. 

5. Conclusions 

We have performed an analysis of the emissivity variations along 
the Galactic plane using the newly released Herschel/Hi-GAL 
data (160 jum < /I < 500 fim) combined with the IRIS 100 fim 
data, at 4' angular resolution. Changes in the emissivity spectra 
are interpreted in terms of the TLS model, which includes the 
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Ttls 


(K) 


Ic (nm) 


A 


Ca 


Full Hi-GAL data 


18.6±1.4 


7.5+5.7 


11.0+9.8 


475* 


Peripheral Hi-GAL data (35° < |/| < 70°) 


18.0+1.2 


7.9±6.5 


18.2+10.8 


475* 


Central Hi-GAL data (|/| < 35°) 


19.4±1.2 


7.2±4.7 


8.3+5.2 


475* 


Solar neighborhood 


17.54 


:0.02 


23.1±22.7 


9.4+1.4 


242+123 










4.8" 


475* 


Compact sources 






5.1+0.1 


3.9+0.1 


1333+68 










10.9" 


475* 


Solar neighborhood + Compact sources 


17.34 


:0.02 


13.4+1.5 


5.8+0.1 


475+20 



Table 1. Parameters of the TLS model and their 1-cr standard deviation; median values of the best-fit parameters for the Galactic 
plane (this work) and best-fit parame ters for the other env ironments (from lParadis et all 1201 lal) . 
* set to the standard value defined in iParadis et al.l (1201 lah . 
**A normalized to ^ = 475. 



DCD process characterized by a coiTelation length (Z^), and the 
TLS effects, whose intensity (A) with respect to the DCD effect 
is left as free parameter Our results can be summarized as fol- 
lows: 

- A 500 fim emissivity excess, with respect to the predictions 
of a modified black-body model with j6=2, has been found in 
the peripheral parts of the Galactic plane (35° < |/| < 70°) 
covered by the data. This excess can represent up to 16% to 
20% of the total emission in some HII regions. 

- The dust temperature, derived from the TLS model, appears 
to be slightly warmer in the central (by 1.4 K, for |Z| < 35°) 
than the peripheral Galactic regions covered by the data. 
Regions near the Galactic center have temperatures in the 
range 17-25 K, whereas the median temperature across the 
rest of the Galactic plane covered is close to 18.6 K. 

- We have found a flattening of the emissivity spectra in 
the range 100-500 //m with increasing dust temperature. A 
model using a mixture of temperatures along the LOS has 
been used to verify wether this could be responsible for the 
observed behaviour. The results strongly suggest that the 
changes in the observed emissivity spectra with dust tem- 
perature cannot be accounted for by an LOS effect only, and 
that they must be caused instead by intrinsic variations in the 
dust properties with environment. These results are indeed 
compatible with the predictions of the TLS model with dust 
temperature, and suggest variations in the degree of amor- 
phisation of the grains, i. e. the disorder at the atomic and/or 
nanometer scales. 

- The 500 fim emissivity excess can be explained by the in- 
tensity produced by the TLS processes. Indeed, the spatial 
variations in the A parameter (i.e. the amplitude of the TLS 
effects) follow the distribution of the excess. 

- Dust properties along the Galactic plane seem to differ from 
those of the solar neighborhood, the excess being smaller in 
the latter than expected from an extrapolation of the Galactic 
plane behavior 

- The spatial distribution of 1^ (i.e. the correlation length of 
the DCD effect) does not present significant variations along 
the Galactic plane. Statistically, is three times lower than 
in the solar neighborhood, although it is close to the value 
obtained for compact sources. 

- Results in the framework of the TLS model indicate that dust 
grains are characterized by a degree of amorphization that 
is larger along the Galactic plane than in the solar neigh- 
borhood. In particular, the degree of amorphization tends to 
increase in the peripheral part of the plane covered by the 
Hi-GAL data. This work shows that specific variations in the 
mechanical structure of the material constituting the grains 



are likely to vary from the Galactic center towards the pe- 
ripheral regions of the plane. 
- Variations in the dust temperature and the 500 fim emissivity 
excess as a function of Galactic longitude appear to correlate 
with the locations of Galactic spiral arms. In the context of 
the TLS model, this may suggest that changes in the mechan- 
ical structure of the grains occur in the spiral arms. 

The combination of Planck and Herschel data will provide 
the opportunity to improve the current constraints on the dust 
properties. 
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Fig. 10. Best-fit model parameters 4- as a function of A, derived for the TLS model, across the entire Galactic plane (left panel), 
in the central (middle panel) and the peripheral regions of the Galactic plane (right panel). The overplotted contours represent the 
density of pixels, with levels at 100, 300, 800, and 1200 for all the Galactic plane and the central regions, and levels at 50, 100, and 
300 for the peripheral regions of the Galactic plane. 
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